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Marek Szydlowski 

Astronomical Observatory, Jagiellonian University, Orla 171, 30-244 Krakow, Poland and 
Complex Systems Research Centre, Jagiellonian University, Reymonta 4, 30-059 Krakow, Poland 

Aleksandra Kurek 

Astronomical Observatory, Jagiellonian University, Orla 171, 30-244 Krakow, Poland 

We develop a simple method of dark energy reconstruction using a geometrical form of the 
luminosity-distance relation. In this method the FRW dynamical system with dark energy is re- 
constructed instead of the equation of state parameter. We give several examples which illustrate 
the usefulness of our method in fitting the redshift transition from the decelerating to accelerating 
. phase as the value of the Hubble function at the transition. 
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^ . I. INTRODUCTION 

> . 

^ ' Recent observations of distant SNIa and other astronomical observations of WMAP, SDSS, Chandra X-ray, etc 
indicate that current Universe is experiencing an accelerated expansion. There are different cosmological models 
which explain this phenomenon by postulating the hypothesis that the Universe is filled by fluid which violates the 
I strong energy condition (so called dark energy) or by postulating that the Friedmann equations do not describe well 
the evolution of the Universe. Due to the lack of a satisfactory cosmological model (while the simplest ACDM model 
seems to be the best one we have still some problems with it) the model-independent reconstruction of the expansion 



I history of the Universe and the nature of dark energy seems to be important. These reconstruction method can 
O be divided into parametric and non-parametric approach Here we consider the reconstruction of dynamics in 
^ , terms of single potential function using the relation between the luminosity distance and the Hubble function. This 
^ ' reconstruction allow us to find the moment of transition from decelerating to accelerating phase of expansion as well 
as to probe the evolution of the dark energy. 



o 



II. RECONSTRUCTION OF THE POTENTIAL FUNCTION OF FRW DYNAMICAL SYSTEM 



> 

. Our Universe in large scale structure is homogeneous and isotropic, therefore its geometry in a good approximation 
' can be described by the Robertson- Walker metric. Assuming general relativity we can characterize evolution of the 
Universe in terms of only single function - scale factor a{t) as a function of cosmic time t. The scale factor obeys the 
Raychaudhuri equation (also called acceleration equation) 



a 1 , ^ 

- = -^P + 3p, 1 
a 



, where dot denotes differentiation with respect the cosmological time t and p — p{t), p = p(t) are energy density and 
^ ■ pressure of the perfect fluid - source of gravity. 

' We assume that both p and p depends on t through the scale factor, i.e. p = p{a{t)) and p = p{a{t)). Then we can 
simply find first integral of Raychaudhuri equation, called as Friedmann equation 

^ , ^ a? p k 

where H is Hubble function, energy density p satisfies conservation condition p — —3H{p-\-p), k — 0, ±1 is curvature 
constant. 

Equation ([T]) can be rewritten to the form of a Newtonian dynamical system 

dV(a) b? . 

where V — V{a) is the potential function which characterizes a model and dot means here differentiation with respect 
rescalled time variable t: t ^ t: \Ho\dt = dr. 
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In this interpretation the Friedmann first integral plays the role of energy integral. For example for the concordance 
ACDM model potential function assumes the following form 

V{a) = -\ J2 ^^,oa~'-'-', (4) 

i—m.k,A 

where Qi o are density parameters for each fluid Qi — defined at the present epoch labelled by index; Wm ~ 

for dust matter, Wk — ~\ for curvature fluid and w — —1 for the cosmological constant; pi = Wi{a)pi in general and 
ao = 1 is chosen for the value of the scale factor at present. 

Let us note that in a general case if equation of state is parameterized by the scale factor a then all dynamics is 
characterized in terms of a single potential function of a (or redshift z: 1 + z — a~^). Of course this function has the 
general form 

Via) = -^n,sia)a^, (5) 
where VlcS is effective energy density parameter for total fluids. 

Since V — Via) ~ _(hJJ^)_s_ ^Yiq potential function Viz) can be immediately reconstructed from SNIa data. Such 
a reconstruction is possible due to a simple relation between Hiz) and luminosity distance for the flat Universe 



d f dLiz) 



dz \ 1 + z 



Hiz). (6) 



Due to the fact that we know ^^(z) for the discrete values of z a kind of smoothing procedure is needed to obtain 
function Viz). This can be done by approximating the luminosity distance function by some fitting formula. Here 

we consider different ansatz for -^^^ = fiO,z) function 

1. polynomial models z) = X]i=o ^i-^^'j where k — 1, ... ,4, 



i=Q' 

2. Fade-type ansatz /(^,z) = where (n,m) = {(1, 1), (1, 2), (1, 3), (2, 1), (2, 2), (3, 1)}, as wefl as fi0,z) 

where x = z + 1 which was proposed in Q , 



x — ai^/x—l~\-ai 



a2+a3\/x+2 — ai 

3. model analysed in fid, z) = 77(1) — rjiy), where 77(2/) = 2a [y~^ + [3y~^ + ^y~^ + 5y~'^ + a] with y = 

Additionally we consider following constraints 1. ^^(z = 0) = 0, 2. iJ(z = 0) = Hq. 

Finally we have chosen this model which is the best one in the light of SNIa data [1, 0, @| in the Baycsian framework 
meaning, i.e. with the greatest value of the posterior probability 0, Q. Here we used the BIG quantity Q as an 
approximation to the evidence and assumed equal values of prior probabilities for all models. 

After such an analysis we conclude that the best parametrization (model) is fiS, z) = 02^^ + z and we have used 
this one in the later analysis. Result of reconstruction the potential function is presented in Figure [1] 

The moment of transition (ot, zt) is related with the maximum of the potential function so it can be easily found. 
Results for the reconstructed potential are as follows OT = 0.59 izT = 0.69), I/(zt) = -0.35 and = 1.98. 

These outcomes can be compared with the values obtained for the potential of the best fit ACDM model ot = 0.57 
izT = 0.75), Vizr) = -0.35 and = 1.47. 

In the next section we consider the scalar field (j) with the potential Uicf)) and show a simple method of obtaining 
Ui(f>ia)) directly from Via). 



III. POTENTIAL OF A SCALAR FIELD FROM THE POTENTIAL OF THE FRW DYNAMICAL 

SYSTEM 

Let us consider a single scalar field with potential t/(0) in the FRW model. They are described by the natural 
lagrangian function 

>C0 = ^5m-^9''</'5''</' -£/(</-). (7) 
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FIG. 1: V{z) for f{6,z) = a2Z^ + z (solid, black line) together with Icr confidence level (dashed lines) with a2 
and for ACDM model (red line) with r2,„,o = 0.27. 
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We assume that a scalar field is homogeneous and minimally coupled to gravity, so kinetic energy term is of the form 

^ 2^ ' — df ^ 

The equation of motion for the field satisfies the Klein-Gordon equation 



(8) 



The previous equation is equivalent to a perfect fluid with energy density and pressure as follows 



(9) 



If we assume that both scalar field as well as potential depends on time through the scale factor then coefficient of 
state = w^{a) can be reconstructed directly from namely 



T + U 
T -U' 



(10) 



And vice- versa if we parametrize ^^(z) by redshift then both kinetic term and potential can be reconstructed from 
simple relation 



U{<P{a)) = i (l-w^{a))p^{a). 
The evolution of quintessence scalar field is governed by adiabatic conservation condition 

/90 + 3H(p^ +P4>) = 0, 

which solution is 

pa / \ " 

-3(l-(-tB^(a)) 



pfiG "exp 



^3 r ^da 



(11) 
(12) 

(13) 
(14) 



where w^{a) 



f (a)d(ln a) 
Jdilna) ■ 



Let us consider flat FRW cosmological model filled by dust matter and dark energy in the form of quintessence 
scalar field. Therefore we have 



Peff = Pm.oa ^ + p^fiE{a), 



(15) 
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where E{a) = a'^ exp [-3 /° ^^da . 

Becausse we assume that matter and dark energy sectors does not interact the potential function is additive, i.e. 



where Vi„(a) = rJm.oi , ^ 



m,0 



V{a) = V^{a) + V4,{a) 
■; V0(a) = fl4,fia'^E{a); Q^^ 



(16) 

It is a simple consequence of the fact that 



Pcl>a 



1 



1 
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= --ii^a . 

From the assumption of fc = we have the Friedmann first integral in the form 

Pm + P0 = Pcff = 3i/^ 
and therefore relation (llip determine dependence of 0(a): 

'c?(?!)\ ^ _ P0(a)(l + i«0(a)) 



da 



[aH{a)Y 



(17) 



(18) 



(19) 



One can check that ([T]) implies that the dark energy coefficient equation of state, ^^(a), is related to the FRW 
potential V as follows 



\ d{\n\V^\) 
d(lna) 



d[\n\{V -V^)\) 



d{\xi a) 



(20) 



where V = —\Q.eff{a)a?^ V,n = —^flm.oo. ^ (see Figure[2]). Hence after substitution ((20l) into (|19p we obtain relation 
(j) = (j){a). 

Let us note that appeared in p9p can be integrated in an exact form in term of the potential V, namely = 



,,oa ^^i"-)^ where 



- 3w0(a) 
(see Figure [2]) which leads to 



ln[|y(a) + ^n,n.,oa-^\] H\V(a - 1) + ifi^.ol) 



ln( 



In ( 



P<i>[ 



P4>.0 



/ V{a) + ^a„,oQ-i 
Vy(a = l) + ir!„^o 



(21) 



(22) 



(see Figure [3]). The analogous procedure of the reconstruction one can perform for the potential of the scalar field 
[/((/)). For this aim we substitute Wcj,{a) from formula ([20|) and P0(a) from ([22|) into (fT2|) . Finally we obtain the 
following formula for reconstruction potential of scalar field as a function of scale factor: 



U{a) = U{a = 1)- 



4F(a = l) + |f7„,o + ^|a. 



(23) 



(see Figure |3]). Both formulas can be rewritten to the form containing redshift z instead the scale factor. 



IV. CONCLUSION 



Different cosmological models can be observationally differentiate in terms of the potential function of the dynamical 
system. We reconstructed the potential function of the FRW dynamical system directly from SNIa observations using 
the relation between the luminosity distance and the Hubble function and the parametric method of smoothing the 
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FIG. 2: w^{z) and W4,{z) for f{0,z) — a2Z^ + z (solid, black lines) together with la confidence level (dashed lines) with 
a2 = -0.2ltoSll, ^rn.o = 0.27 and w{z) for the ACDM model (red lines). 




FIG. 3: ^ and J^^l^^-^ for /(^, z) = 022^^ + 2: (solid, black lines) together with la confidence level (dashed lines) with 
aa = -0.2ltoMh ^rn.o = 0.27 and for ACDM model (red line). 

function dL^z). As one concludes the shape of the reconstructed potential is similar to the shape of the ACDM 
model potential for small redshift [z < 1). We compute the moment of transition using the reconstructed potential as 
well as the potential of the ACDM model. We also consider the cosmological model with a single, homogeneous and 
minimally coupled to gravity scalar field with the potential U (0). We derive the relation between 'w^{a) / p^{a) /U{a) 
function and the potential function of the FRW dynamical system V{a) which allow us to reconstruct the evolution of 
such functions without making assumptions on the form of the and p^. We compare evolution of such quantities 
with the evolution predicted by the ACDM model. 
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